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Taking into account a neutral massive scalar field minimally coupled to gravity, 
in a Robertson- Walker metric, it is shown that when the final state is connected with 
the initial one by means of a Bogoliubov transformation, which does not include 
the single-mode rotation operator, the mean value of created particles is conserved. 
When the rotation operator is considered, it is still possible to use the approach of 
single-mode squeezed operators and get the entropy as the logarithm of the created 
particles. 
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In the work by Gasperini and Giovannini [1] it is shown that the large amount of 
entropy of our present universe can be explained by the mechanism of cosmological 
pair creation from the vacuum, as a consecuence a background of squeezed cosmic 
gravitons is produced. In our case we will study the contribution of a neutral 
massive scalar field to the entropy, without any graviton present, in a classical 
curved background. In order to do that we will use the formalism introduced by 
Parker [2] which let us define the general time dependent Bogoliubov transformation 
in a Robertson- Walker metric. In that reference a consistent method to obtain the 
annihilation-creation operators at each time is given. In this work we will show 
that the general time depending Bogoliubov transformation, introduced in ref. [2], 
transforms the "in" vacuum to a squeezed state. 

Ref. [1] is based on the approach of ref. [3] where a linearized graviton theory is 
used. In that case the graviton modes satisfy the Klein-Gordon equation correspond- 
ing to a massless scalar field minimally coupled to the gravitational background. In 
the present case we use a massive scalar field. In both cases the procedure is similar. 
In ref. [4] we can see that still in the massless case the particle creation number is 
different from zero, due to the kind of coupling (with the conformal coupling the 
particle creation is zero [5]). 

Let us now rewrite the expresions of ref. [2] that we need in the following. The 
minimally coupled scalar field satisfies the field equation 

(V M ^ + m 2 V = (1) 

(with n = 0,1,2,3). 

For an unbounded space-time, as we use, a convenient technical method is the 
discretization in the momentum space, which is implemented as in ref. [2] by means 
of the boundary condition </?(x + nL) = (p(pc,t), where n is a vector with integer 
Cartesian components and L a longitude which goes to infinity at the end of the 
calculation (in our case it is not necessary because the results used do not depend 
on L). Then a general solution of eq. (1) can be expanded in the form 

V?(x, t) = 5}a k (*)0k(a;) + a\{tyt>\{x)] (2) 

k 

where ^ 

6k(x) = — ^expi(kx- / W(k, t')dt') 

V y ' (La(t)) 3 / 2 v / 2W F 1 Jto 

whith W a real function of k = |k| and t. The operators a^(t) and dk'(i) are 
time dependent and satisfy the commutation relations 
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M*), ak'(t)] = 0, ^kW.a+k'C*)] = 0, [a k (£),aV(£)] = 5 k)k , (3) 
The operators act on the vacuum |0, t > in the form 

a k (*)|0,*>=0, a t k (t)|0,t>= |l k ,t > 
Moreover we define the operators 

^4 k : = o k (t = ti) ujzt/i i > ti > to 
V :=a k t (t = ti) 

Also we define the vacuum |0>:=|0,t = ti>on which the operators defined before 
act. These operators are related with the time dependent ones by the following 
Bogoliubov transformation 

a k (£) = e l7 " (k ' t) coshr(k,t)A k + e l7 ' 3(k ' t) sinhr(k,t)A t _ k (4a) 

a_ k f (t) =e" l7 ' 3(k ' t) sinhr(k,t)A k + e- l7Q(k ' t) coshr(k,t)A t _ k (46) 
where the functions 7 Q and r satisfy the obvious conditions 

7a (k,t 1 ) = 0. 
r(Mi) = 0. 

Replacing eqs.(4) in eq. (2) and all of that in eq. (1) we obtain the following 
system of equations (see ref. [6]) 

(1 + tanhrcosr)M + 2W r 7 a = (5a) 

Msinr + 2WY = (56) 
with r := 7a + 7/3 - 2 f* Wdt' and 

M := -l(W/W) + \{W/Wf - ^{a/af - ha/a)+u 2 - W 2 (5c) 

£i 4r 4r 
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where uj 2 = k 2 /a 2 + m 2 , "a" is the scale factor of the universe and the dot is 
the cosmic time derivative. 

Following ref. [2] we can rewrite eqs. (4) in the form 



a k = e F A k e- F (6) 

with the convenient notation 



-k 

k= — oo 



F k ,_ k = ±r(k, t)[e- 2iX ^A k A_± - e^^A^A^^] + z 7a (k, t){A\, A k } 

where 2A(M) = 7a (k, t) + 7/3 (k, t). 
Defining also F as: 

Using that F commutes with -F_ k k , ^k,-k and A k , we can rewrite eq.(6) in 
the form 

a k = exp(F_ kjk + F kj _ k )A k exp(-F_ k)k - F kj _ k ) (7) 

Let us introduce the definitions of two-modes squeezed and single-mode rota- 
tion operators as in ref. [7]. Using the commutation relation between A k and A]J 
and taking into account the isotropy of the space by means of the assumption 

r(-k) = r(k) 
la (-k) = 7a (k) 
A(-k) = A(k) 

we can rewrite eq.(7) in the form 

a k (t) =S{ r -, \)R^{-2 la )sC-, A) J R k (-2 7a )A k 

S\ r -, \)tf^(-2 la )S\ r -, A) J Rt k (_2 7a ) ( 8 ) 



where 



S( r -,X) = exp{^[e- 2 ^^)A k A_ k - e 2 ^A\A^ k }} 



is the two-mode squeezed operator, and 

#k(-27a) = exp{2z 7a (k,t)A t k A k } 



is the single-mode rotation operator. 
From ref.[7] we can use the relation 

R(2 1 )S(r/2,X + 1 ) = S(r, A)i?(2 7 ) 

i?(-2 7 )A kJ Rt(_2 7 )=e- 2 ^ k 

Therefore we obtain 

a k (£) = e" 2 ^(r/2, X)S(r/2, A + 7a )A k ^(r/2, A + 2 7a )S t (r/2, A + 7a ) (9) 

In eq. (9) we can see that the deviation from the pure squeezing transformation 
is given by the function 7q . 

We will study now the particular case 

7 a = (10) 



then we have 

a k (t) = S(r,X)A k S\r,X) (11) 

with S(r, A) a two-mode squeezed operator. Following ref.[l] we can now ex- 
pand eq.(ll) in terms of only single-mode squeezed operators. In order to do that 
we can introduce the bosonic operators 6 k and 6_ k by means of the definition 

-4k = ;-^(& S i 5 (k)k - isig(k)b_ stg(k)k ) (12) 

where sig(k) is the sign function, defined by 

sig(-k) = -sig(k) 

and by a set of directions which are choosen as positive. We can separate the states 
into the following disjoint sets {k} = II k U Il_ k , where 

II k = {k so that k.u > when k.u ^ 0; 

k.v > when k.v ^ and k.u = 0; k.(v x u) > when k.u = k.v = 0} 
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with u and v two arbitrary mutually perpendicular directions in the k-space. 
Moreover 

Il_ k = {k so that k = -k', with k' e II k } 

Now we can define sig(\s) = 1 if k e Ilk and -1 if k e Il_k. In what follows we 
suppose that k G Ilk. 

Then we can rewrite the two-mode squeezed operator as the product of single- 
mode operators: 

S(r,\) = S 1+ (r,\)S 1 -(r,\) (13) 

with 



5 1± =exp[^(e- 2a 6 ± k 2 -e 2lA 6t ±k 2 )] 

Replacing eqs (12) and (13) in eq. (11) we obtain 

a k = -±=(S 1+ b k S\+ - iSt-b-kS*!-) (14) 
V2 

Using now the procedure of ref.[l] we can factorize, the wave function ty zk in 
the space (x+,x-) which corresponds respectively to Ilk and Il_k, in the form 

^ zk (x+x-) =< x+x-\z k >= tf 1 zk {x + )^ 1 zk (x-) 

It is not necessary to repeat the procedure of ref. [1], which finally permits us 
calculate the reduced density operator by means of 



Pk = J dx + dx-\x + X- >< x + X-\z k >< z k \x + X- >< X + X-\ 
and the entropy using 

S = —Trplnp 

Therefore the entropy associated to the k-mode can be calculated (as is shown 
in ref.[l]), by means of the approximate equation (valid when rik >> 1): 

S(k)~lnn k (15) 

But we will see that this expression gives us only the entropy for the "in-out" 
process. We have not a dynamic equation for S(k), because the assumption given 
by eq.(10) is valid when is a constant. 
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We will analize now the physical consequences of the no rotation condition 
given by eq. (10). If 7 Q is the the null function, eqs (5) turn to 



(1 + tanhrcosT)M = 



(16a) 



MsinT + 2Wf = 



(166) 



We have two options in order to satisfy eqs. (16): 

a) One possibility is by means of the following condition on the function M (see 
eq. 5c), for any time t>t±: 



from eq. (16b) it turns out that the squeezing parameter must be a constant 
in time. Therefore the mean value of created particles between the times t\ and 
t is independent of time, because from eqs (4) we have 



Then, particles only can be created if in equal number are annihilated. 

The condition (17) is really a differential equation for the function W, which 
can be considered as the "frequency" of the normal modes, in an expansion of the 
field as the one given by eq. (2). 

It is interesting to note that, when the process of interaction between the 
matter field and the geometry is considered as "in-out", with ti = ti n and t = t out , 
the conditions on W(t out ) and W(t out ), that in the literature [8] are known as the 
diagonalization of the hamiltonian, satisfy automaticaly eq.(17). So a process where 
the final result is a hamiltonian with the functional form of independent oscillators, 
produces squeezed states. 

b) The other possibility to satisfy eq. (16), is the requirement 



M = 



(17) 



?T-k =< 0|a T kak|0 >= sinh rk 



(18) 



cos r = — coth r 



(19) 



Now performing the derivative of eq. (19) and replacing in eq. (16) we obtain 



Putting eq.(20) in terms of the cothr and replacing in eq. (19), we get the 
following integro-differential equation for the function W and the phase 7^: 



M[W] 1/2 cos[j(3 -2 Wdt'] + {2W(^f3-2W) + M[W}} 1/2 = (21) 

J t 



From eqs (20) and (21) we can get a relation between the average number of 
particles "n" and the function W. In order to do that we can rewrite eq. (20) in 
the form 

* = + 2W 

Integrating we have 

ft 



7/3 



\L^ At '^L wit '^ (22) 



with \x an integration constant. 

Replacing now eq.(22) in eq.(21) and using the assumption that M 7^ 0, in 
order to make the division by M 1 / 2 , we get: 



1 r* M 1 1/2 



Since the module of the cosine must be greater than one and n is a possitive 
integer, eq.(23) only makes sense when n = 00. Therefore the only non trivial 
possibility to solve eqs (16) is the one given in a). In that case it is possible to 
calculate the entropy by means of eq. (15), as in ref. [1]. But, when the dynamics 
of the particle creation given by eq.(5) is taken into account, the transformation of 
the operators is given by eq. (9). 

We will analize now, if it is still possible to apply the same approach to calculate 
the entropy, in the case when the dynamics, given by eqs (5), is not altered, 

The transformation between operators given by eq. (9) is not a pure squeezing 
one. However from eq.(9) we can obtain, using the unitarity of the operators; 

flk(*)kk>=0 (24) 

with 

W >:= S(r/2, A + 7«)S(r/2, A + 2 7a )|0 > (25) 
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then \z > is a squeezed state, because it results of the application of two 
successive squeezing transformation. 

As in ref. [1] we can introduce the wave function ^ Z ^(x + ,X-) :=< x + ,X-\z >, 
which, from eq. (24), satisfies the equation 

a k *,k = (26) 



In order to make explicit eq. (26), we can take into account eq. (4a), replacing 
eq. (12) there, then we have the operator a^(t) as a function of 6k, 6-k, and its 
hermitic conjugate: 

v / 2a k (t) = e lla coshr6 k + e %1<3 sinhr6 t k - i(e lla coshr6_ k + e %l13 sinhr6 t _k) (27) 
We can now introduce the following condition in eq.(27) 

7/3 = 7« (28) 



It was used that 7/3 does not appear either in eq. (5) or in the mean value of 
particle creation, therefore it is a "gauge condition" that we can choose in order to 
obtain a Gaussian distribution for the wave function in full analogy with ref. [1]. 

Taking into account the commutation relations [6, b'] = 1 and [x,p] = i. Calling 
to the superfluctuant operator related with the operator 6k and X- to the one 
corresponding to 6_k, we can write 

b^ = ^=-(x + + d x+ ) , b- k = ^=-(x-+d x _) (29) 

and analogously for the conjugated ones. Replacing these expressions in eq.(27) 
it results 

i[coshr(x + + d x+ ) — sinhr(x + — d x+ )]^ + [coshr(a;_ +d x _) — sinhr(a;_ — d x _)]^ = 



Introducing the factorization 

^(x + x-) := ^ 1 (x + )^ 1 (x-) 



For both variables the function ty 1 satisfies the equation 

coshr(x + dx)^ 1 = sinhr(x — dx)^ 1 (30) 
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Then ^ is a Gaussian distribution and all the calculation performed in ref. 
[1] is available, so the entropy can be approximated by the logarithm of the mean 
value of the created particles. Therefore in a way similar to that in ref. [9] it is 
possible to apply this formalism to calculate the cosmological entropy production. 

It is interesting to note that we can get, for a large particle creation number, 
the same dependence with the entropy, if we use the expression of the entropy 
corresponding to a bosonic gas in flat space time (see ref. [10]), but with the mean 
value of particle creation coming from the dynamics of the curved space time. 
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